Data-constrained reionization and its effects on cosmological parameters 
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We perform an analysis of the recent WMAP7 data considering physically motivated and viable 
reionization scenarios with the aim of assessing their effects on cosmological parameter determina- 
tions. The main novelties are: (i) the combination of CMB data with astrophysical results from 
quasar absorption line experiments; (ii) the joint variation of both the cosmological and astrophys- 
ical [governing the evolution of the free electron fraction Xf,{z) ] parameters. Including a realistic, 
data-constrained reionization history in the analysis induces appreciable changes in the cosmological 
parameter values deduced through a standard WMAP7 analysis. Particularly noteworthy are the 
variations in Vlih^ = 0.02258lo:ooo56 (WMAP7) vs. Q.bh^ = 0.02183 ± 0.00054 (WMAP7 + AS- 
TRO), and the new constraints for the scalar spectral index, for which WMAP7 + ASTRO excludes 
the Harrison- Zel'dovich value ris = 1 at > 3a. Finally, the e.s. optical depth value is considerably 
decreased with respect to the standard WMAP7, i.e. Te — 0.080 ± 0.012. We conclude that the 
inclusion of astrophysical datasets, allowing to robustly constrain the reionization history, in the ex- 
traction procedure of cosmological parameters leads to relatively important differences in the final 
determination of their values. 

PACS numbers: 



I. INTRODUCTION 

It is well known from a large set of astrophysical ob- 
servables that after primordial recombination (which oc- 
curred at a redshift of z ^ 1100) the universe "reionized" 
at a redshift z > 6. It is common practice in Cosmic 
Microwave Background (CMB) studies to parametrize 
the reionization as an instantaneous process occurring at 
some redshift Zr, with 4 < z,. < 32, and to marginalize 
over Zr when deriving constraints on the other cosmo- 
logical parameters. In the absence of any precise astro- 
physical model of the reionization process, the electron 
ionization fraction Xe{z) is parametrized by Zr in the fol- 
lowing way: Xe{z) = 1 for z ^ z^ (possibly Xe{z) — 1.08 
or Xe{z) — 1.16 for z < 3 in order to take into account the 
first and second Helium ionization) and Xe{z) < 2 x 10~^ 
for z > Zr in order to join the ionization fraction value 
after the recombination. In the following we will refer 
to this parametrization as "sudden" or "instantaneous" 
reionization. With this choice of parametrization there 
exist a one-to-one relation between the redshift of sud- 
den reionization Zr and the electron scattering optical 
depth Tg. The most recent constraints on the optical 
depth that come from the analysis of the Wilkinson Mi- 
crowave Anisotropy Probe team on their seven-year data 
(WMAP7), in which it is assumed a sudden reionization 
scenario, is Tej = 0.088 ± 0.015. However, as already no- 
ticed, e.g. in [ll, and further emphasized by our previous 
works ( 2] and 3] ) , the assumption of a general reioniza- 
tion scenario could affect the extraction of the constraints 
of cosmological parameters. In particular, we studied the 
effects of non-instantaneous reionization on the two prin- 
cipal inflationary parameters (the scalar spectral index 



of primordial perturbations and the tensor-to-scalar 
ratio parameter r), and on the optical depth t^. The 
method used in the above cited works to describe a gen- 
eral reionization scenario, developed in Ref. [l|, is based 
on a principal components (PC) analysis of the reion- 
ization history, Xe{z) . PCs provide a complete basis 
for describing the effects of reionization on large-scale E- 
mode polarization spectrum. Following Ref. [l|, one can 
treat Xe (z) as a free function of redshift by decomposing 
it into its principal components: 



Xe{z) = xl{z) + ^ m^S^{z), 



(1) 



where the principal components S^{z) are the eigenfunc- 
tions of the Fisher matrix describing the dependence 
of the polarization spectra on Xe{z)', the are the 
PC amplitudes for a particular reionization history, and 
xl{z) is the WMAP fiducial model for which the Fisher 
matrix is computed and from which the PCs are ob- 
tained. Therefore the amplitude of eigenmode /i for 
a perturbation around the fiducial reionization history 

SXe{z) = Xe{z) — xl{z) IS 
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dz Sfj_{z)SXe{z) 



(2) 



In Refs. and [1] we made use of the publicly avail- 
able S'^(z) functions and varied the amplitudes for 
the first five eigenfunctions (i.e. for /i = 1,...,5). The 
principal components were computed only in the range 
of redshifts z S [6 — 30] . 

In what follows we refer to this parametrization of 
reionization as the "Principal Components" (PC) reion- 
ization. Since the ionization fraction is bounded in 
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< Xe{z) < 1 (neglecting helium reionization and the 
small residual ionized fraction after recombination) in 
the range of redshifts in which PCs are defined, it is nec- 
essary to impose some limits on the amplitudes of the 
eigenmodes of equation ([2]) to let the reionization frac- 
tion be within these limits, if only for the definition of 
reionization fraction. In Ref. ^j] the authors find the 
ranges of values for the amplitudes compatible with 
Xk{z) G jo, 1] for all the redshifts in range of interest. In 
[31 and Q we performed a Monte Carlo Markov Chains 
analysis assuming a flat prior on (only) the ranges of val- 
ues of the amplitudes whose linear combination with 
the function give a Xe{z) in the allowed range. These 
values are reported in left part of Table HV] and are la- 
beled "PC Bounds". 

However, these limits for the values of the PC ampli- 
tudes are a necessary but not sufhcient condition for the 
reionization fraction to lie in < Xe{z) < 1. In fact, as 
noticed also by ^J, if any rrifj, violates those bounds Xe{z) 
is guaranteed to be unphysical in some redshift range, but 
the opposite is not true, because the full reionization his- 
tory depends on the linear combinations of the product 
of the amplitudes times their corresponding PC principal 
component. Indeed, even if all the amplitudes satisfy 
the bounds reported in Table ITVl Xe{z) could assume an 
unphysical value for some redshifts. To overcome this 
potential problem, we have added in the version of the 
cosmomc package used in Q and Q the condition that 
the value of Xe{z) computed at each step of a Markov 
Chain must be in the range < Xe{z) < 1 for every z. In 
these studies, this was the only "physicality" condition 
imposed on the possible reionization history. However, 
experimental data gathered in the last few years can be 
used to discard at least some of the possible Xe{z) histo- 
ries on well understood (astro)physical grounds. It is now 
possible to use reionization histories that are physically 
motivated and tested with known probes of the reioniza- 
tion epoch, such as the Gunn-Peterson optical depth, or 
the distribution in redshift of the Lya emitters. 

In this work we adopt the results of a well-tested semi- 
analytical reionization model proposed in Refs. Q and 
(in what follows we will refer to this model as the CF 
model). This model takes into account a large number 
of parameters and physical processes that are involved in 
modeling reionization, including (e.g.) the radiative and 
chemical feedbacks of the first sources of ionizing light 
on the evolution of the intergalactic medium (IGM), and 
constrain the model by comparing it with a variety of 
observational data, such as the redshift evolution of Ly- 
man Limit Systems (LLS), the IGM temperature and the 
cosmic star formation density. Thus we will be able to 
build up an ensemble of reionization histories that is more 
robust from both the theoretical and the observational 
point of view, rather then rely on purely phenomenologi- 
cal, albeit model- independent, parameterization schemes 
as the PCs. 

We will combine the CF model with a standard ACDM 
cosmological model and we perform a Monte Carlo 



Markov Chains analysis of the joint CMB and reioniza- 
tion data. We will thus be able to test the impact of 
considering a detailed physical model for reionization on 
the constraints of the cosmological parameters, and con- 
versely to test the dependence of the CF model on the 
underlying cosmological model. 

At the end of such analysis we will moreover derive the 
subsequent constraints on the amplitudes of the reion- 
ization principal components (applying directly the 
equation [J) . By construction then, these limits on the 
values of amplitudes of the principal components will be 
compatible and constrained both by the CMB and by the 
astrophysical probes of the reionization process. 

The main objectives of the present work are then: 

• Verify the impact of considering a data-constrained 
and realistic reionization model on the determina- 
tion cosmological parameters. 

• Verify the impact on the constraints of the reion- 
ization parameters produced by variations of the 
cosmological parameters, i.e. refraining from fixing 
them a priori from the most updated best fit values 
of the WMAP experiment. 

• Obtain the PC amplitudes from the allowed 
reionization histories. 

As such an analysis with combined cosmological pa- 
rameters characterizing the background evolution of 
the universe and astrophysical parameters modeling the 
reionization history has not yet been made, it is worth- 
while to explore their mutual implications on the extrac- 
tion of the constraints of the two ensemble of parameters. 

II. ANALYSIS 

The details of the CF model are summarized in Ref. 
8]; in the present work we assume the following settings: 

• We consider here a flat ACDM cosmology described 
by a set of cosmological parameters: 

{flm,^bh'^,h,as,ns}, (3) 

where rim is the total matter density relative to the 
critical density, ftbh^ is the baryonic matter density, 
h is the reduced Hubble parameter Ho — lOOh, ag is 
the r.m.s. density fluctuation in spheres of radius 
8h^^ Mpc and Ug is the scalar spectral index of 
primordial perturbations. We want to stress that 
these cosmological parameters are considered here 
as free parameters, so that they are not assumed a 
priori, as in Q. 

• The CF reionization model contains additional 
three free parameters. These are enjn = 
[e*/osc]ii,iii, the product of the star- forming effi- 
ciency (fraction of baryons within collapsed haloes 
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Parameter 


Mean 


95% c.l. limits 




0.2733 


[0.2260, 0.3305] 




0.2184 


[0.0208, 0.0229] 


h 


0.6984 


[0.6553, 0.7422] 


Us 


0.9579 


[0.9330, 0.9838] 


CT8 


0.7941 


[0.7434, 0.8491] 


en 


0.0037 


[0.0016, 0.0067] 


£iii 


0.0165 


[0.0000, 0.0398] 


Ao 


3.0152 


[1.0000, 5.1739] 


Te 


0.0803 


[0.0625, 0.1042] 


Zr 


6.7469 


[5.8563, 8.2000] 



TABLE I: Mean and 95% c.l. constraints on the cosmological, astrophysical and derived parameters obtained with the reion- 
ization parametrized with the CF model of reionization. 



going into stars) and the fraction of photons es- 
caping into the IGM /osc for PopII and PopIII stars; 
the normahzation Aq of the ionizing photons mean 
free path (see Ref. |8j] for details). In what follows 
we refer to these three parameters as the "astro- 
physical" parameters, to distinguish them from the 
five "cosmological" ones described described in the 
previous point. 

• The ranges of variation adopted for the three free 
astrophysical parameters are en G [0;0.02], eni G 
[0;0.1], Ao e [1;10]. 

• The observational data used to compute the like- 
lihood analysis are (i) the photo-ionization rates 
Fpi obtained using Lya forest Gunn-Peterson op- 
tical depth observations and a large set of hydro- 
dynamical simulations @ and (ii) the redshift dis- 
tribution of LLS diVLL/dz in the redshift range of 
0.36 < 2; < 6 The data points are obtained 
using a large sample of QSO spectra. For details, 
see Ref. @. 

• In order to make the analysis self-consistent, the 
WMAP7 constraint on the total electron scatter- 
ing optical depth is not considered in this analy- 
sis. This prevents a possible loophole in our analy- 
sis: WMAP7 constraints on have been obtained 
using the assumption of instantaneous reionization 
sX z = Zr- Once this idealized evolution of Xe{z) 
is dropped (this paper), the value of Tg must be 
a byproduct of the new analysis rather than being 
inserted artificially as an external constraint into 
it. Moreover, as already pointed out in Ref. Q, 
the CMB polarization spectra are sensitive to the 
shape of the reionization history and considering 
a more general reionization scenario could lead to 
a tighter optical depth constraint than derived by 
WMAP7 [2j. 

• Finally, we impose the prior that reionization 
should be completed by z = 5.8 to match the flux 
data of Lya and Ly/3 forest. 



With these hypotheses we have then modified the 
Boltzmann CAMB code [l3| to incorporate the CF model 
and performed a MCMC analysis based on an adapted 
version of the public available MCMC package COSMOMC 



|12| . Our basic data set is the seven-yr WMAP data [11 1 
(temperature and polarization), on top of which we add 
two "astrophysical" datasets, i.e. the LLS redshift evo- 
lution, dN^L/dz Ref. 7], and the Gunn-Peterson optical 
depth measurements presented in Ref. [g} . To extract the 
constraints on free parameters from such combined data 
set we consider a total likelihood function L cx exp(— £) 
made up by two parts: 



C 



T 



obs 



T 



th 



(4) 



where £' refers to the WMAP7 likelihood function and 
is computed using the routine supplied by the WMAP 
team; Ja represents the set of A^obs observational points 
referring to Gunn-Peterson optical depth LLS distribu- 
tion data; finally, are the corresponding observational 
error-bars. We constrain the free parameters by maxi- 
mizing L with fiat priors on the allowed parameter ranges 
and the aforementioned prior on the end of reionization 
at z = 5.8. 

The Monte Carlo-Markov Chain convergence diagnos- 
tics are done on 4 chains applying the Gelman and Rubin 
"variance of chain mean" /"mean of chain variances" R 
statistic for each parameter. We considered the chains to 
be converged at i? — 1 < 0.03. 



III. RESULTS 

The results of the MCMC analysis described above are 
summarized in Table HI where we list the marginalized 
posterior probabilities at 95% confidence level (c.l.) er- 
rors on the free cosmological and astrophysical param- 
eters. We also report the constraints for two derived 
parameters: the electron scattering optical depth and 
reionization redshift Zr. to be intended as the redshift 
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Parameter 


WMAP7 


WMAP7 + PC 


WMAP7 + ASTRO 


h 

0-8 


0.266 ±0.029 

0.710 ±0.025 
0.963 ±0.014 
0.801 ± 0.030 


0.243 ± 0.032 
0.02321 ± 0.00076 
0.735 ± 0.033 
0.994 ± 0.023 


0.273 ±0.027 
0.02183 ± 0.00054 
0.698 ± 0.023 
0.958 ±0.013 
0.794 ± 0.027 


Te 
Zr* 


0.088 ±0.015 
10.5 ± 1.2 


0.093 ± 0.010 


0.080 ± 0.012 
6.7 ±0.6 



TABLE II: Comparison of the 68% c.l. posterior probabilty contraints obtained for different parametrizations of reionization. 
* The Zr parameter has a different definition in the different reionization scenarios (see text for details). 



Parameter 


WMAP7 + ASTRO Mean 


WMAP7 ± ASTRO 95% c.l. limits 


CF Mean 


CF 95% c.l. limits 


£// 


0.0037 


[0.0016, 0.0067] 


0.003 


[0.001, 0.005] 




0.0165 


[0.0000, 0.0398] 


0.020 


[0.0000, 0.043] 


Ao 


3.0152 


[1.0000, 5.1739] 


5.310 


[2.317, 9.474] 


Tel 


0.0803 


[0.0625, 0.1042] 


= 0.088 ±0.015 


= 0.088 ± 0.015 


Zr 


6.7469 


[5.8563, 8.2000] 


6.762 


[5.800, 7.819] 



TABLE III: Comparison between the mean value and the 95% c.l posterior constraints between the present work (WMAP7 + 
ASTRO) and the CF model, Ref. H (MCF). 



at which the reionization is 99% complete. In Table 
im we show the 68% c.l. constraints obtained by the 
WMAP team for the standard 6-parameter ACDM model 
("WMAP7") and the constraints obtained on the cosmo- 
logical parameters from the present analysis ("WMAP 7 
+ ASTRO"). 

As we can see from the Table HI] the results of our work 
mildly differ from the WMAP 7 results for the parameters 
of the standard ACDM model. The most sensitive pa- 
rameter for the presence of the "astrophysical" datasets 
(LLS and Gunn-Peterson data) is ni,h^ whose mean val- 
ues in the two cases differ by more than a standard devi- 
ation from each other. It is important to note that even 
when considering a complex reionization history imply- 
ing three new parameters the errors remain practically 
the same as in the standard case. 

Table |TT] reports the results obtained in Q for the 
WMAP7 dataset with the PC reionization ("WMAP7 + 
PC"). This method produces two main differences with 
respect to the WMAP7 + ASTRO case: the first is re- 
lated to the constraints obtained for Us- In Q the con- 
straints for the scalar spectral index were compatible 
with Tis = 1, i.e. the Harrison-Zel'dovich (HZ) primor- 
dial power spectrum, when instead WMAP7+ASTR0 
excludes the value = 1 at > Scr. The second differ- 
ence concerns Te in the two cases: for WMAP 7 + PC 
this quantity is in the range Te = 0.093 ± 0.010, while 
the WMAP+ASTRO case gives a mean value lower by 
> l-fj, i.e. Te = 0.080±0.012. Note that in the WMAP7 
+ PC case we did not consider constraints on the as pa- 
rameter, so in Table HTl the corresponding value is missing. 

There is a caveat in comparing the constraints ob- 
tained on Zr- Indeed, in the WMAP7 case Zr is the red- 
shift at which the universe undergoes an instantaneous 



and complete reionization process. In the more realistic, 
extended reionization scenarios considered here instead, 
Zr is defined as the redshift at which the IGM is 99% 
re-ionized by volume. With this clarification in mind, 
WMAP7+ASTR0 results predict 5.8 < Zr < 8.2 at 95% 
c.l. (see Table m. 



In Table Hill we report the 95% c.l. posterior probabil- 
ity constraints for the reionization parameters e//, em 
and Ao obtained in the present work (WMAP7 + AS- 
TRO case, cosmological parameters free to vary) com- 
pared to those obtained in Ref. ^ in which the cosmo- 
logical parameters were fixed to the WMAP7 best fit val- 
ues (CF case). Figure [T] shows the comparison between 
the best-fit model for the Xe{z) evolution for the two 
cases of WMAP7 + ASTRO and CF. For the WMAP7 
+ ASTRO case, full hydrogen reionization is not only 
achieved earlier than in the CF model, but the evolu- 
tion is faster, resulting in an initially lower Xe{z) above 
z — 8. These differences are entirely induced by the fact 
that we have now allowed the cosmological parameters 
to vary together with the astrophysical ones, but they 
are relatively small. The fact that the astrophysical pa- 
rameters do not show much dependence on cosmology 
is understandable because the cosmological parameters 
affect the reionization process mostly through structure 
formation. The next obvious step is to include large scale 
structure information in the analysis. In conclusion, in- 
cluding astrophysical datasets in the analysis seems to 
lead to relatively important effects on the extraction of 
the cosmological parameters. 
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FIG. 1: Ionization histories for the best-fit model for the 
two cases WMAP7+ ASTRO (red dotted solid curve) and CF 
(green solid curve) @]. 



A. PC amplitude reconstruction 

For each reionization history allowed by the MCMC 
likelihood analysis, we use eq. ^ to reconstruct the 
amplitudes of the first five PC amplitudes, m^, with 
fi = 1...5. By construction now, the amplitudes not 
only fulfill the necessary physicality conditions (see Sec. 
1) but also they are compatible with the additional astro- 
physical data sets considered in this analysis, i.e. the Lya 
Gunn-Peterson test and the LLS redshift distribution. 

In Fig. [2] we show the two dimensional 68% and 99% 
c.l constraints for the amplitudes obtained here com- 
pared with those obtained in [2[ for which we show the 
two dimensional 68% and 95% c.l. distributions for each 
of the cases considered. We choose to report the 99% c.l. 
instead of the usual 95% c.l. limits to be as conservative 
as possible in showing the reionization histories allowed 
by the MCMC likelihood analysis. The color (layer) code 
is the following: in pink (top layer) there is the case 
WMAP7 + ASTRO considered in the present work. In 
the background there are the cases considered in in 
blue is the WMAP7 case (bottom layer) , in red (next 
layer up) is the case called "CMB All" ( i.e. WMAP7 
+ ACBAR -I- BICEP-K QUAD + BOOMERanG), green 
(next layer) is CMB All + LRG-7 and yellow (next layer) 
is simulated Planck data. 2] considered an ensable of 
CMB dataset along with WMAP7, and also we forecasted 
future constraints from the Planck experiment, simulat- 
ing a set of mock data with a fiducial model given by the 
best fit WMAP5 model with the following experimental 
noise: 




Nf 



w 



-1/2 



^K-rad 



exp 



£{£ + l){0 



FWHM 



/rad)2 



iln2 



(5) 



FIG. 2: 68% and 99% reconstructed c.l. constraints 
for the values of the PC amplitudes computed from CF 
model and eq. [5] (top layer, pink). Background con- 
tours refer to 68% and 95% c.l. constraints obtained in 
01 with the PC reionization for WMAP7 (bottom layer, 
blue), WMAP7+QUAD-HACBAR-HBICEP (CMBAU, next 
layer up, red), CMB All -I- LRG-7 (next layer, green) and 
simulated Planck data (next layer, yellow), respectively. 



where w"^/^ is the temperature noise level (a factor \/2 
larger for polarization noise) and 9 is the beam size. For 
the Planck mission we use w^^^ = 58^K and ^fwhm = 
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Parameter PC Bounds Astrophysical Bounds 



mi 


[-0.1236,0.7003] 


[-0.1229, -0.0866] 


7712 


[-0.6165,0.2689] 


[-0.2594, 0.0002] 


ms 


[-0.3713, 0.5179] 


[0.0763, 0.2941] 


7714 


[-0.4729, 0.3817] 


[-0.2107, -0.1080] 


7715 


[-0.3854, 0.4257] 


[0.0418, 0.1319] 



TABLE IV: Ranges of variation for tiie amplitudes of the principal component, in the case of the Principal Components and 
in the case of the 99% c.l. reconstructed amplitudes of the present analysis (see text for details). 



7.1' equivalent to expected sensitivity of the 143 GHz 
channel. 

The region spanned by PC amplitude values is much 
smaller than that allowed by when the PC bounds only 
are imposed. The 99% c.l. constraints values are re- 
ported in the right part of the Table pVI) ( "Astrophysical 
Bounds" ) . As seen from Table (|IV|) the amplitudes of all 
the principal components (except for 7712) obtained with 
the above procedure are constrained at 99% c.l. to take 
a definite sign, negative for mi and 7714 and positive for 
7773 and m^. Moreover, even if the 99% c.l. upper bound 
of 7712 is positive, this second amplitude is mostly con- 
strained to be always negative. These results are in qual- 
itative agreement with [^j who also found that the same 
amplitude signature, albeit with errors large enough that 
the 95% c.l. bounds encompass values of both possibile 
signs. 

IV. CONCLUSIONS 

With the aim of constraining the evolution of cosmic 
reionization, we have extended previous work based on 
the use of Principal Components analysis. The main nov- 
elty of the present work is represented on one hand by 
complementing available CMB data with additional as- 
trophysical results from quasar absorption line experi- 
ments, as the Gunn-Peterson test and the redshift evo- 
lution of Lyman Limit Systems. In addition, we have 
for the first time explored the effects of a joint variation 
of both the cosmological {ilm, ri&ft-^, h, cts, Tig) and astro- 
physical (eii, em, Ao, see Sec HI] for their physical mean- 
ing) parameters. Note that, differently from the vastly 



used approach in the literature, we do not impose a pri- 
ori any bound on the electron scattering optical depth Tg, 
which instead we calculate a posteriori. This is to prevent 
a possible loophole in the calculation, as the WMAP de- 
termination of such quantity is based on the assumption 
of an instantaneous reionization which we do not make 
here. 

Including a realistic (i.e physically motivated) reion- 
ization history in the analysis induces mild changes in 
the cosmological parameter values deduced through a 
standard WMAP7 analysis. Particularly noteworthy 
are the variations in flbh^ = 0.02258t^;|5J^^^^ (WMAP7) 
vs. Qbh^ = 0.02183 ± 0.00054 (WMAP7 + ASTRO), 
and the new constraints for the scalar spectral in- 
dex, for which WMAP7-I- ASTRO excludes the Harrison- 
Zel'dovich value ti^ = 1 at > 3(t. Finally, the e.s. op- 
tical depth values is considerably decreased with respect 
to the standard WMAP7, i.e. = 0.080 ± 0.012. We 
conclude that inclusion of astrophysical datasets, allow- 
ing to robustly constrain the reionization history, in the 
extraction procedure of cosmological parameters leads to 
relatively important differences in the final determination 
of their values. 
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